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CHAPTER  I 

INTRODUCTION 

Aircraft  designed  for  flight  at  low  supersonic  Mach 
numbers  have#  for  a long  time,  employed  the  swept  wing  with 
subsonic  leading  edge  to  prevent  the  appearance  of  shock 
waves  in  the  flow  field  (at  least  ahead  of  the  trailing  edge) 
and  thereby  avoid  the  losses  associated  with  them.  At  higher 
Mach  numbers,  say  eibove  3,  such  a wing  would  be  so  highly 
swept  that  its  low  speed  performance  is  unacceptable.  One 
must  therefore  go  to  either  a variable  geometry  wing  or  a 
wing  with  a supersonic  leading  edge.  In  the  latter  case,  it 
is  desirable  from  the  standpoint  of  lifting  efficiency  to 
have  the  shock  waves  attached  to  the  leading  edge  and  con- 
fined to  the  flow  over  the  lower  surface.  Wings  supporting 
a flow  of  this  type  are  known  collectively  as  waveriders  (1) . 

The  general  problem  of  calculating  the  flow  field  about 
an  arbitrary  three-dimensional  body  in  a supersonic  stream 
is  unmanageable.  Thus  another  approach  is  usually  taken  with 
waveriders,  the  inverse  method,  wherein  one  starts  with  a 
simple  known  flow  field  and  constructs  a body  to  fit  it  by 
replacing  stream  surfaces  with  the  solid  walls  of  the  body. 
The  flows  which  have  been  used  in  this  context  are  two- 
dimensional  wedge  flow,  two-dimensional  Prandtl-Meyer  expan- 
sion and  euci-symmetric  cone  flow. 
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The  first  published  report  of  a body  of  this  type  is  by 
G.  I.  Maikapar  (2)  of  the  U.S.S.R.  in  which  a pyramid  whose 
base  is  a polygon  with  concave  (reentrant)  sides  is  surrounded 
by  n regions  of  wedge  flow  with  the  plane  oblique  shocks 
spanning  the  n corners  of  the  polygon.  A lifting  body  de- 
rived from  wedge  flow  was  first  reported  by  T.  Nonweiler  (3,4) 
and  is  actually  1/n^^  of  the  Maikapar  body.  This  body,  the 
caret  wing,  has  an  inverted  V cross  section  from  which  it 
received  its  name  (Figure  1) . J.  W.  Flower  (5)  has  described 
upper  expansion  surfaces  derived  from  a Prandtl-Meyer  expan- 
sion and  the  combination  of  these  with  caret  wing  lower 
surfaces.  A survey  of  much  of  the  work  done  in  Great  Britain 
using  the  inverse  method  has  been  presented  by  J.  Seddon  (1) 
and  A.  Spence. 

All  bodies  derived  from  a known  flow  field  necessarily 
have  a "design  point",  the  free  streeun  Mach  number  of  the  de- 
riving flow  and  the  angle  of  incidence  of  the  body  to  that 
flow.  At  other  flight  conditions  the  simple  two-dimensional 
or  euci-symmetric  flow  will  not  be  maintained  and  some  approxi- 
mate method  must  be  used  to  describe  the  off -design  behavior. 
J.  Venn  (6)  and  J.  W,  Flower  have  given  a qualitative 
description  of  the  shock  patterns  for  off-design  caret  wings 
and  have  shown  that  exact  solutions  other  than  at  the  design 
point  exist  for  certain  of  these  wings.  The  flow  field  of  a 
body  similar  to  the  reentrant  pyramid  but  with  blunt  rather 


than  sharp  corners  and  a shock  attached  only  at  the  apex  has 
been  solved  numerically  by  B.  Ziph  (7)  using  a finite 
difference  approximation  to  the  equations  governing  conical 
flow.  The  approximate  method  used  in  this  study  will  be 
described  in  Chapter  II. 

This  study  is  concerned  only  with  the  caret  wing  wave- 
rider.  The  lower  surface  of  the  caret  wing  can  be  described 
by  the  three  parameters  6,  0 and  $ (Figure  1).  Delta  is  the 
angle  between  the  internal  rib  and  the  free  stream  direction. 
Theta  is  the  angle  between  the  plane  spanning  the  leading 
edges  and  the  free  stream  direction.  Phi  is  the  angle  between 
the  two  free  stream  surfaces  containing  the  leading  edges; 
i.e.,  the  emgle  between  the  projection  of  the  leading  edges 
onto  a plane  perpendicular  to  the  free  stream.  The  volume  in 
6,  0,  $ space  of  all  possible  caret  wing  lower  surfaces  is 
shown  in  Figure  2 as  the  cylinder  with  triangular  ends  ABC 
and  DEF.  The  ranges  of  the  parameters  are  0°  £ <I>  £ 180®, 

0°  £ 0 £ 90®  and  0®  £ 6 £ 0.  Note  that  in  the  limit  $ -►  180® 
the  caret  wing  becomes  a two-dimensional  wedge. 

Caret  wings  for  which  a design  Mach  number,  exists 

form  a more  restricted  class.  At  the  design  point  a plane 
oblique  shock  spans  the  leading  edges  and  5 and  0 are  the 
stream  deflection  angle  and  shock  angle,  respectively,  of  the 
corresponding  wedge  flow.  For  any  given  0 there  is  a maucimum 
6,  occurring  at  infinite  Mach  number,  above  which  wedge  flow 
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does  not  have  an  attached  shock  solution.  Thus  M,  will 

des 

exist  only  for  those  values  of  6 less  than  or  equal  this 
maximum.  In  Figure  2,  the  volume  of  all  caret  wings  that 
have  a design  point  is  the  cylinder  with  end  surfaces  AGB 
and  DHE . 

This  volume  can  itself  be  divided  into  two  parts,  de- 
pending on  whether  the  wedge  flow  corresponding  to  any 
6,  0 pair  is  a weak  or  strong  shock  solution.  For  wedge 
flow,  the  strong  shock  is  not  observed  in  practice  and  it  is 
therefore  not  expected  to  occur  in  the  caret  wing  flow  field 
either.  Only  caret  wings  derivable  from  wedge  flow  of  the 
weak  type  will  be  considered  in  the  r^ainder  of  this  paper. 

In  Figure  2,  the  volume  in  5,  0,  ♦ space  of  these  caret 

wings  is  the  cylinder  with  end  surfaces  AGIB  and  DHJE.  The  J 

' 

terms  "weak"  and  "strong"  are  used  above  in  the  conventional 

sense  and  should  not  be  confused  with  their  use  later  on  in  1 

discussions  about  the  component  of  the  free  stream  velocity 

normal  to  the  leading  edge. 

The  caret  wing  body  studied  here  will  be  completed  with 
an  upper  expansion  surface  defined  in  a manner  analogous  to 
the  lower  surface  and  a base  surface  normal  to  the  free  stream 
direction  (Figure  3) . Just  as  the  lower  surface  is  defined 
by  the  parameters  6,  0 and  $,  the  upper  surface  will  be  de- 
fined by  6',  0 and  <I>  where  6'  is  the  angle  between  the  upper 
surface  rib  and  the  free  stream  direction  and  0 and  $ are  the 
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seune  as  for  the  lower  surface.  The  coimnon  values  of  0 and  4* 
guarantee  that  the  leading  edges  of  the  lower  and  upper  sur- 
faces coincide.  To  reduce  the  number  of  independent  variables, 
' will  always  be  taken  equal  to  3°  less  than  6 . As  a re- 
sult, the  thickness  ratio  will  be  about  0.05. 

The  main  concern  of  this  study  is  the  estimation  of  the 
lift  coefficient,  drag  coefficient  and  lift-to-drag  ratio  of 
caret  wings  and  the  identification  of  their  dependence  on 
body  geanetry,  Mach  number  and  Reynolds  number.  Mach  numbers 

and  Reynolds  numbers  considered  range  from  1.5  to  5.0  and  from 
5 8 

10  to  10  respectively.  The  aerodynamic  coefficients  are 
calculated  by  combining  estimates  of  the  upper,  lower  and 
base  surface  pressures  and  estimates  of  the  skin  friction  on 
the  upper  and  lower  surfaces.  Calculations  of  the  pressures 
on  the  upper  and  lower  surfaces  are  based  on  supersonic,  in- 
viscid,  perfect  fluid  gas  dynamics  and  calculations  of  base 
pressures  rely  on  published  experimental  data.  The  skin 
friction  calculations  are  for  flat  plates  with  zero  pressure 
gradient. 


A compeurison  of  the  performance  of  caret  wings  to  delta 
wings  will  also  be  made. 


CHAPTER  II 


A MODEL  OF  THE  CARET  WING  FLOW  FIELD 

2.1  Analysis  of  the  Caret  Wing 
Flow  Field  at  Design  Conditions 

It  is  well  known  that  the  velocity  component  parallel  to 
an  oblique  shock  wave  passes  through  unchanged  and  that  the 
shock  acts  as  a normal  shock  to  the  normal  velocity  component. 
Except  for  relations  involving  ratios  of  static-to-total 
conditions,  the  changes  in  properties  across  an  oblique  shock 
are  the  same  as  for  the  normal  velocity  components  and  a 
normal  shock.  On  this  basis  the  use  of  normal  shock  tables 
is  extended  to  oblique  shock  waves  (8) . More  generally,  one 
can  take  any  oblique  shock  flow,  add  or  subtract  a velocity 
component  of  any  magnitude  which  is  tangent  to  the  shock  and 
calculate  - based  on  the  new  velocities,  shock  angle  and 
stre2uti  deflection  angle  - the  change  in  properties  across  the 
shock  in  the  original  flow. 

An  interesting  phenomenon  associated  with  this  calcula- 
tion procedure  is  that,  while  the  changes  in  properties  across 
the  shocks  are  the  same  for  both  flows,  the  shock  type  ("weak" 
or  "strong")  may  change.  In  a simple  case,  for  example,  to 
find  the  pressure  ratio  across  a weak  oblique  shock  in  two- 
dimensional  wedge  flow  one  might  calculate  the  Mach  nximber  of 
the  velocity  component  normal  to  the  shock  and  look  up  the 
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pressure  ratio  in  a normal  shock  table.  However,  the  normal 
shock  is  of  the  strong  type  (the  weak  solution  corresponding 
to  the  normal  shock  is  a Mach  wave)  and  thus  the  question 
arises,  whether  to  use  the  weak  or  strong  solution.  Clearly, 
the  spatial  location  of  the  shock  cannot  change.  In  the 
example  under  discussion,  the  Mach  number  component  used  in 
the  table  look  up  is  normal  to  the  shock  and  the  solution 
which  preserves  this  orientation,  the  strong  solution,  is  the 
desired  one. 

A slightly  different  application  of  the  concept  of  sub- 
tracting velocity  components  tangent  to  a shock  wave  is  found 
in  the  analysis  of  the  swept  wedge  of  infinite  span.  The 
component  of  the  free  stream  velocity  which  is  parallel  to 
the  leading  edge  of  the  wedge  will  also  be  tangent  to  any 
shock  attached  to  the  leading  edge.  Therefore,  this  velocity 
component  can  be  subtracted  from  the  flow  field  by  a Galilean 
transformation  of  coordinates  leaving  an  upstream  velocity 
which  is  normal  to  the  leading  edge  and  which  is  deflected  by 
the  wedge  through  the  angle  the  wedge  makes  as  measured  in  a 
plane  perpendicular  to  the  leading  edge.  This  is  just  two- 
dimensional  wedge  flow.  If  the  normal  velocity  component  is 
supersonic  and  the  wedge  angle  normal  to  the  leading  edge  is 
small  enough,  weak  and  strong  attached  shock  solutions  are 
mathematically  possible.  In  practice  the  weak  shock  is 
always  observed  in  two-dimensional  wedge  flow  and  the 
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existence  of  crossflow  in  the  swept  wedge  problem  will  not 
alter  this,  even  though  the  strong  shock  with  respect  to  the 
normal  velocity  component  is,  in  some  cases,  weak  with  re- 
spect to  the  free  stream.  If  the  normal  velocity  component 
is  subsonic  or  the  wedge  angle  is  too  large,  no  attached 
shock  solution  exists. 

These  ideas  will  now  be  applied  to  the  caret  wing  flow 

field.  The  caret  wing's  lower  surface  is  formed  from 

sections  of  two  swept  wedges  which  are  joined  at  the  vertical 

symmetry  plane  of  the  wing.  At  the  design  point,  a plane 

shock  spans  the  leading  edges  and  deflects  the  free  stream 

through  the  angle  6,  the  angle  between  the  free  streaim  and 

the  internal  rib.  In  the  same  manner  as  the  swept  wedge  of 

infinite  span,  the  component  of  the  free  stream  velocity 

parallel  to  either  of  the  caret  wing's  leading  edges, 

must  be  tangent  to  the  design  point  shock  and,  therefore, 

can  be  subtracted  from  the  free  stream,  V . The  shock  can 

00 

then  be  analyzed  using  the  component  of  the  free  stream 
velocity  normal  to  the  leading  edge,  and  the  wedge  angle 
normal  to  the  leading  edge,  Pj^,  just  as  well  as  from  the 
point  of  view  of  and  6. 

Consider  a family  of  caret  wings,  all  derivable  from  the 


1 

'.1 


1 


seune  two-dimensional  wedge  flow.  They  will  have  the  angles  6 
and  0 in  common  while  4>  varies  from  zero  to  180  degrees 
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(Figure  4) , A single  design  Mach  number  determined  by 
<5  and  Q will  also  be  common  to  the  members  of  this  family. 

In  the  limit  ♦ * 180°  the  caret  wing  is  a two- 
dimensional  wedge.  The  velocity  component  normal  to  the 
leading  edge  is  equal  to  the  free  stream  velocity  and  the 
wedge  angle  normal  to  the  leading  edge  is  equal  to  6.  As 
mentioned  in  the  introduction,  only  caret  wings  whose  design 
point  shocks  are  weak  with  respect  to  the  free  stream  are 
being  considered.  So  when  * 180°  the  design  point  shock 
is  also  weak  with  respect  to  V^. 

In  the  limit  4 = 0°  the  velocity  component  noinnal  to 
the  leading  edge  is  also  perpendicular  to  the  design  point 
shock.  Thus  the  shock  must  be  strong  with  respect  to  V|_. 

As  explained  earlier,  this  shock  is  never  found  on  the  swept 
wedge  of  infinite  span  and  therefore  doubt  is  cast  on  the 
existence  of  the  caret  wing's  design  point  flow  as  4 approaches 
zero.  However,  the  fact  that  the  strong  shock  is  never  ob- 
served on  an  isolated  swept  wedge  is  not  sufficient  grounds 
for  dismissing  the  possibility  of  its  occurrence  here,  since 
each  leading  edge  of  the  caret  wing  lies  within  the  Mach  cone 
of  influence  of  the  opposite  wing  panel.  More  will  be  said 
about  the  tenability  of  these  strong  shock  solutions  in  the 
next  section. 


Evidently,  at  some  value  of  4>  between  zero  and  180  de- 
grees the  design  point  shock's  type  with  respect  to  ^ under- 
goes a transition  from  weak  to  strong.  The  value  of  $ at 
which  this  occurs,  ,must  be  characterized  by  a smooth 

transition  between  types  because  physically  the  shock  is  the 
seime  regardless  of  *,  only  the  velocity  conponent  being 
analyzed  changes  with  $.  For  any  given  Mach  number  there  is 
a maucimum  angle,  through  which  a streamline  can  be 

deflected  by  a single  oblique  shock  wave  and  at  which  the 
weak  and  strong  shocks  coincide  (the  governing  cubic  equation 

has  a double  root) . Thus  is  that  value  of  $ where 

trans 

r^ax  • This  is  illustrated  in  Figure  5 for  the 

particular  caret  wing  family  6 = 5.0°  and  0 = 34.0°,  where 


Hneix'  ' 


and  P -r  are  plotted  as  functions  of 

max  -1- 


^max  -i’j.  ' "trans*  ''2.29  . 

Caret  wings  whose  design  point  shocks  are  weak  with  re- 
spect to  can  now  be  further  classified  according  to  whether 
the  shock  is  weak  or  strong  with  respect  to  Wings  with 

♦ greater  th£m  -have  the  weak  solution  and  those  with  ♦ 

^rans 

less  than  strong  one.  * trans  plotted  a?  a func- 

tion of  0 with  5 as  a parameter  in  Figure  6,  Using  the  curve 
corresponding  to  the  value  of  5 one  is  interested  in,  caret 
wings  lying  to  the  left  of  the  curve  have  design  point  shocks 
that  are  weak  with  respect  to  and  those  lying  to  the  right. 
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strong  shocks.  The  range  of  9 for  a given  5 varies  from  a 

minimum  when  M,  = ® to  a maximum  when  P (M,  ) * 6. 

des  max  des 

It  is  interesting  to  plot,  for  the  flow  normal  to  the 

leading  edge  of  a family  of  caret  wings,  the  pressure  ratio 

of  both  the  weak  and  strong  shocks  as  a function  of  4>.  This 

is  done  for  the  family  5 » 5.0°  and  0 » 34.0°  at  M = ^(jes 

♦ 2.02  in  Figure  7.  Starting  at  4>  = 180°  and  down  to 

♦ = '^trans*'  42.29°  the  weak  shock  gives  the  design  point 
pressure  ratio  which  is  a constant  independent  of  4>.  At 

♦ « '*’trans^^®  strong  shock  picks  up  where  the  weak  shock 

leaves  off  and  continues  the  design  point  pressure  ratio 
down  to  = 0°.  The  weak  and  strong  solutions,  patched  to- 
gether at  ifr  « give  the  flow  field  the  family  of 

caret  wings  was  derived  from. 

2.2  Extension  of  the  Swept  Wedge 
Analysis  to  Off-Design  Conditions 

Schlieren  photographs  taken  in  wind  tunnel  studies  of 
Maikapar  bodies  at  Mach  nuiiibers  below  indicate  that  a 

convex  curved  shock  attached  to  the  leading  edges  replaces  ij 

the  plane  shock  found  at  the  design  point  (9) . For 

and  Venn  (6)  and  Flower  suggest  a con- 

vex shock  composed  of  a curved  central  section  joined  to  plane 
shocks  emanating  from  the  leading  edges  (Figure  8) . The  shock 
is  plane  outside  the  Mach  cone  from  the  caret  wing's  apex  and 
is  the  weak  shock  calculated  using  and  . 
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For  > Mjjgg/  Venn  (6)  and  Flower  suggest  either  a 
bifurcated  shock  pattern  or  a pattern  where  plane  shocks 
from  the  leading  edges  intersect  at  the  vertical  symmetry 
plane  and  are  followed  by  a complex  system  of  shocks  further 
inside  the  wing  (Figure  8) . In  either  case  the  leading  edge 
shocks  are  plane  outside  the  Mach  cone  from  the  wing’s  apex, 
are  inclined  inside  the  plane  spanning  the  leading  edges  and 
are  the  weak  shock  calculated  using  and  . 

These  complicated  shock  patterns  will  be  approximated 
here  by  two  plane  shocks,  one  from  each  leading  edge,  that 
intersect  at  the  wing's  vertical  symmetry  plane  (Figure  8). 

It  was  shown  in  the  previous  section  that  the  on-design 
caret  wing  can  be  treated  as  a swept  wedge.  The  model  of  the 
caret  wing's  off-design  behavior  adopted  here  will  be  to 
extend  the  use  of  the  swept  wedge  calculations  to  all  condi- 
tions, except  for  $ < Within  this  approximation,  the 

pressure  on  the  caret  wing's  lower  surface  is  a constant  de- 
termined by  the  pressure  ratio  across  the  leading  edge  shocks. 

The  model  just  described  for  the  caret  wing's  lower 
surface  flow  field  gives  the  exact  solution  under  the 
assumptions  of  steady,  supersonic,  inviscid,  perfect  fluid 
aerodynamics  when  the  free  stream  Mach  number,  M^,  is  equal 
to  the  design  Mach  number.  It  also  gives  the  exact  solution 
in  the  limits  'll  = 180°  or  6 = 0°. 
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Under  the  same  assumptions,  the  model  gives  a shock 
tangent  to  the  actual  shock  at  the  leading  edges  regardless 
of  the  influence  of  the  opposite  wing  panel.  This  is  be- 
cause the  Mach  number  cind  stream  deflection  angle  at  the 
leading  edges  are  known  and  are  sufficient  to  determine  the 
shock.  Furthermore,  the  shocks  given  by  the  model  will  also 
be  correct  outside  of  the  Mach  cone  from  the  apex  of  the 
caret  wing.  It  should  be  noted  that  this  model  of  the  flow 
field  violates  the  boundary  condition  that  there  be  no  flow 
across  the  vertical  symmetry  plane. 

The  average  pressure  on  the  lower  surface  of  a caret 
wing  is  approximated  by  the  pressure  downstream  of  the 
model's  leading  edge  shocks.  Useful  bounds  on  the  magnitude 
of  the  error  in  this  approximation  are  not  known.  However, 
something  can  be  said  about  the  sign  of  the  error.  A method 
of  determining  the  magnitude  of  the  error  will  be  suggested 
in  Chapter  V. 

For  cases  where  M < M , , the  model's  shocks  are  in- 

flo  aes 

dined  outside  of  the  plane  spanning  the  leading  edges. 
Therefore,  they  direct  the  flow  away  from  the  vertical 
symmetry  plane  and  expansions  will  be  present  in  the  actual 
flow  which  curve  the  shock  toward  the  internal  rib  so  that  it 
is  perpendicular  to  the  symmetry  plane  at  their  intersection 
(6).  Thus  the  model  gives  too  high  a pressure  in  this  case. 
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For  cases  where  > M , , the  model's  shocks  are 

“ des 

inclined  inside  the  plane  spanning  the  leading  edges. 
Therefore  they  direct  the  flow  toward  the  symmetry  plane 
and  further  compression  will  be  present  in  the  actual 
flow  (6).  In  this  case  the  model  gives  a pressure  below 
the  correct  value. 

More  can  now  be  said  about  the  strong  shock  solutions 

for  $ < At  the  design  Mach  number  it  was  possible 

to  make  a smooth  transition  from  the  weak  to  strong  cases. 

This  transition  occurred  at  / where  T is  equal  to 

trans  max 

and  the  weak  and  strong  solutions  coincide.  is  a 

geometric  property  of  the  caret  wing  and  does  not  change 
with  M^.  However,  is  a monotonically  increasing  func- 

tion of  M . If  M is  raised  above  M,^„,  then  there  will  no 
longer  be  any  value  of  $ for  which  F is  equal  to  F,  and 
no  smooth  transition  from  the  weak  to  the  strong  cases  is 

possible.  In  Figure  5 the  entire  curve  for  F vs.  $ is 

nidx 

shifted  upwards  when  M > M , and  F - F never  goes  to 

00  des  mcLX  j-  ^ 

zero.  It  is  not  likely  that  a discontinuous  change  in  shock 
pattern  or  lower  surface  pressure  will  occur  with  an 
infinitesimal  change  in  body  geometry  or  flight  conditions 
and  for  this  reason  the  strong  shock  solutions  must  be  con- 
sidered untenable.  Other  possibilities  include  a weak  shock 
at  the  leading  edge  followed  by  a complicated  system  of 
shocks  further  inside  the  wing,  unsteady  behavior  or  detach- 
ment. Venn  (6)  and  Flower  suggest  the  first  of  these  on  the 


grounds  that  weak  shocks  aure  always  observed  on  two- 
dimensional  wedges.  Because  of  the  uncertainty  here,  caret 
wings  with  design  point  shocks  that  are  strong  with  respect 
to  Vj_  will  not  be  considered  in  calculations  of  lift,  drag 
and  lift-to-drag  ratio. 

One  additional  generalization  about  the  caret  wing's 

off-design  behavior  will  be  made.  When  is  reduced  below 

M, , the  entire  curve  for  r vs.  $ in  Figure  5 is  shifted 

des  max 

downward.  This  leaves  a neighborhood  about  in  which 

r is  less  than  F.  and  consequently,  the  shock  is  detached. 

IllaX  ^ 


2.3  The  Upper  Surface  Flow  Field 

The  model  for  the  upper  surface  flow  field  is  analogous 
to  the  model  for  the  lower  surface.  Corresponding  to  the 
shock  on  the  lower  surface  is  a Prandtl-Meyer  expansion  of 
the  flow  normal  to  the  leading  edge  for  the  upper  surface. 

The  component  of  the  free  stream  velocity  normal  to  the 
leading  edge,  is  the  scime  as  before  and  the  expansion 
angle  is  calculated  in  the  same  way  as  except  that  6'  is 
used  instead  of  6.  The  average  pressure  on  the  upper  surface 
of  a caret  wing  is  approximated  by  the  pressure  downstream  of 
the  Prandtl-Meyer  expansion  at  the  leading  edge.  Since  the 
expansion  at  the  leading  edge  directs  the  flow  toward  the 
vertical  symmetry  plane,  additional  compression  will  take 
place  in  the  actual  flow  and  the  model's  estimate  of  the  upper 
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surface  pressure  must  be  too  low.  This  model  also  violates 
the  boundary  condition  of  no  flow  across  the  symmetry  plane. 


2.4  Estimate  of  Skin  Friction 

An  approximation  of  the  skin  friction  drag  consistent 
with  the  constant  pressure  surfaces  of  the  upper  and  lower 
surface  flow  field  models  is  employed  here.  The  calculation 
is  for  a two-dimensional  flat  plate  with  zero  pressure 
gradient.  The  Mach  numbers  and  Reynolds  numbers  downstream 
of  the  lower  surface  shock  and  upper  surface  expansion  are 
used  as  the  free  stream  conditions  in  the  calculation  of  the 
skin  friction  on  the  lower  and  upper  surfaces  respectively. 

A power  law  for  the  variation  of  viscosity  with  temperature 
valid  in  the  range  300°R  to  900°R  is  used  to  calculate  the 

Reynolds  numbers  downstream  of  the  shock  and  expansion  (8) . ^ 

Skin  friction  estimates  for  both  laminar  and  turbulent 
boundeiry  layers  are  made. 

The  local  value  of  the  skin  friction  coefficient  for 
laminar,  incompressible  flow  based  on  distance  from  the  lead- 
ing edge  in  the  direction  of  the  external  flow  (parallel  to 
the  internal  rib  for  the  lower  surface  and  parallel  to  the 
upper  surface  rib  for  the  upper  surface)  is  integrated  over 
the  lower  and  upper  wetted  surfaces  to  obtain  the  average 
skin  friction  coefficients  on  the  lower  and  upper  surfaces. 
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These  are  multiplied  by  the  ratio  of  compressible  to  in- 
compressible skin  friction  coefficients  based  on  an  adiabatic 
wall,  Prandtl  number  equal  to  one  and  exponent  in  the  vis- 
cosity law  equal  to  0.765  given  by  Schlichting  (10)  for  the 
Mach  number  range  zero  to  ten. 

Schlichting  (10)  found  an  empirical  equation  which,  in 
5 9 

the  range  5 x 10  < Re^j^  < 10  , fits  the  exact  formula  for 

the  turbulent,  incompressible  skin  friction  coefficient 
based  on  a logarithmic  velocity  profile.  This  is  integrated 
over  the  lower  and  upper  wetted  surfaces  to  give  the  average 
skin  friction  coefficients  on  the  lower  and  upper  surfaces. 
These  are  then  multiplied  by  the  ratio  of  compressible  to 
incompressible  skin  friction  coefficients  from  the  theory 
of  R.  E.  Wilson  as  reported  by  Schlichting  (10) . This  ratio 
is  valid  for  the  Mach  number  range  zero  to  ten. 

2.5  Estimate  of  Base  Pressure 

The  calculations  of  caret  wing  base  pressure  rely  on 
published  experimental  data.  The  main  source  of  this  data 
is  Chapman  (11),  Wimbrow  and  Kester's  wind  tunnel  study  of 
blunt  trailing  edge  wings  at  Mach  numbers  1.5,  2.0  and  3.1 
for  both  laminar  and  turbulent  boundary  layers.  Base  pressure 
data  taken  on  the  vertical  fin  of  the  X-15  at  Mach  5 has  been 
used  to  extend  the  Mach  number  range  of  the  turbulent  boundary 
layer  data  (12).  The  work  of  Rom  (13),  Gadd  (14)  and  Korst 
(15)  have  been  compared  to  Chapman's  data.  Base  pressure  data 
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taken  on  Maikapar  bodies  at  the  Aerophysics  Laboratory  (16,17) 
have  also  been  compared  to  the  other  data.  Chapman's  graphs 
of  base  pressure  ratio  vs.  ratio  of  boundary  layer  thickness 
to  trailing  edge  thickness  appecu:  in  Figure  9 with  the  other 
authors'  data  added  for  comparison.  The  Maikapar  body  base 
pressures  tend  uo  be  higher  than  the  other  data.  This  is 
apparently  due  to  the  lower  aspect  ratio. 

Chapmcin  (11)  found  that  base  pressure  data  at  a given 
free  stream  Mach  number  is  correlated  to  the  ratio  of  boundary 
layer  thickness  at  the  trailing  edge  to  trailing  edge  thick- 
ness. This  variable  depends  on  the  thickness  ratio  and 
Reynolds  ntamber  based  on  wing  chord.  The  thickness  ratio  of 
a caret  wing  is  constant  at  all  spanwise  locations.  The 
free  stream  Reynolds  number  based  on  half  the  length  is  used 
as  an  average  value  over  the  span. 

Straight  lines  have  been  fit  to  Chapman's  (11)  turbulent 
boundary  layer  base  pressure  curves  and  a line  through  the 
X-15  data  point  and  parallel  to  Chapman's  M^  = 3.1  line  has 
been  used  for  M^  = 5.0.  Curves  have  been  fitted  to  Chapman's 
laminar  boundary  layer  base  pressure  curves.  At  intermediate 
Mach  numbers  linear  interpolation  between  the  curves  is  used. 


CALCULATION  OF  THE  LIFT  COEFFICIENT,  DRAG 
COEFFICIENT  AND  LIFT-TO~DRAG  RATIO 


I 

3 . 1 Introduction 

There  cure  six  independent  vairiables  in  this  problem. 

Four  of  them,  <S,  6',  0 and  $,  define  the  cauret  wing's 
geanetry  and  incidence  to  the  free  stream.  The  other  two, 
the  free  stream  Mach  number,  M_,  and  the  free  stream 

Reynolds  number,  Re^,  specify  the  flight  conditions.  For  |i 

t. 

convenience,  the  free  stream  Reynolds  number  has  been  based  ii 

I;. 

on  the  square  root  of  the  plan  area,  S,  so  that  by  holding  | 

Re^  constant  wings  of  the  same  plan  area  rather  than  length  | 

can  be  compared.  A simple  trigonometric  formula  which 
relates  the  Reynolds  numbers  based  on  /S  and  length  is  | 

presented  in  Section  3.6.  i 

■j 

J 

The  lift  coefficient,  C^^,  drag  coefficient,  Cj^,  and  j 

j 

lift-to-drag  ratio,  L/D,  are  to  be  calculated  as  functions  j 

of  the  six  independent  variadsles.  Five  cases  are  considered; 

1,  no  skin  friction  or  base  drag;  2,  with  Iciminar  skin  fric- 
tion but  no  base  drag;  3,  with  turbulent  skin  friction  but 
no  base  drag;  4,  with  laminar  skin  friction  and  base  drag; 

5,  with  turbulent  skin  friction  and  base  drag. 


J 
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Given  a set  of  values  for  the  independent  variables,  a 
preliminary  procedure  must  be  undertaken  to  make  sure  the 
caret  wing  belongs  to  the  group  being  considered  in  this 
study.  To  be  a member  of  this  group  must  exist,  the 

design  point  shock  must  be  weak  with  respect  to  the  free 
stream  and  ♦ must  be  greater  than  or  equal  to 

trans 

If  the  above  conditions  are  satisfied,  the  calculation 
of  Cjj  and  L/D  can  proceed.  The  equations  for  inviscid, 

compressible  flow  are  taken  from  NACA  Report  1135  (8) . The 
ratio  of  specific  heats,  c /c  , is  assumed  equal  to  1.4. 

3.2  Mach  Number  and  Wedge  Angle 
Normal  to  the  Leading  Edge 

Formulas  for  Mj_  and  will  be  needed  in  three  phases 

of  the  calculation  of  C_ , C_  and  L/D.  They  are  necessary  in 

the  determination  of  pressure  ratio  across  the 

lower  surface  shock  and  the  pressure  ratio  across  the  upper 
surface  expansion.  In  the  derivation  which  follows,  the 
Mach  number  will  be  treated  as  a vector. 

Let  the  x-eixis  of  a cartesian  coordinate  system  be 
aligned  with  the  free  stream  direction  and  the  y and 
z-aoces  be  vertical  and  horizontal  respectively.  Also,  let 
the  origin  of  coordinates  be  located  at  the  apex  of  a caret 
wing  of  unit  length  (Figure  3) . Point  A is  the  caret  wing 
apex,  point  B is  the  left-hand  wing  tip  and  point  C is  the 
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intersection  of  the  internal  rib  with  the  base.  Vector  B 
is  a vector  pointing  from  the  origin  to  point  B. 

The  (X/Y/z)  coordinates  of  points  A,  B and  C are 


A * 

(0,0,0) 

(3.2.1) 

B = 

(l/b2  ^^3 ) 

(3.2.2) 

C = 

(1,C2»0) 

(3.2.3) 

where 


^2 

= - tan 

0 

(3.2.4) 

^3 

= tan  0 

tan 

2 

(3.2.5) 

c 

0 

= - tan 

6. 

(3.2.6) 

The  ccanponent  of  parallel  to  the  left-hand  leading 
edge,  will  be  calculated  by  forming  the  dot  product  of 

and  a unit  vector  parallel  to  the  leading  edge.  The  free 
stream  Mach  number  in  vector  form  is 

Moo  * M^i  = (^^,,0,0)  . (3.2.7) 

A unit  vector  parallel  to  the  leading  edge  is 


(if  ^2 ^^2  ^ 
/l+b2+b3 


(3.2.8) 
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The  Mach  number  parallel  to  the  leading  edge  is 


M 


cf 


b| 


III 


2 ,2 


(3.2.9) 


l+b2+b3 


and  its  magnitude  is  given  by 


M 


cf 


I "of  I 


V 


l+b2^b3 


(3.2.10) 


Calculation  of  the  Mach  number  normal  to  the  left-hand 


leading  edge  is  now  simply  a matter  of  subtracting  M^^  from 


M 


= ««  - ^Cf 


.2.^2 


2 J (^32'^^3'  "^2'  "^3^  (3.2.11) 

l+b2+b3 


= 


4 2 2 4 2 2 

^2b^i  + bj+b^+b^ 


5 T 

1 + b^  + b^ 


(3.2.12) 
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The  wedge  angle  normal  to  the  left-hand  leading  edge, 
r^,  is  the  angle  between  and  the  plane  containing  the 

left-hand  lower  surface  of  the  caret  wing  (Figure  3) . 

It  will  be  found  using  the  dot  product  of  and  a unit 

vector  nonnal  to  the  lower  surface,  ABCj_. 

The  equation  of  a plane  with  normal  vector  ()l,m,n)  is 


ax  + my  + nz  = c.  (3.2.13) 

Equation  (3,2.13)  can  be  solved  for  a,  m,  n and  c subject 
to  the  condition  that  the  plane  contain  the  points  A,  B and  C 
and  that  vector  (a,m,n)  be  of  unit  length.  The  solution  is 


a 


m 


-C2in 


n-1/2 


C2  + 1 + 


^2"^2 


n 


b 


3 


2 


)m 


c » 0. 


(3.2.14) 

(3.2.15) 

(3.2.16) 

(3.2.17) 


i) 


Thus  a unit  vector  normal  to  the  caret  wing's  left-hand  lower 
surface  is 


2 . , . 

C2  + 1 + — ) 


-1/2 


^2~^2 


(3.2.18) 


ABC,  is  the  inward  normal  due  to  the  selection  of  the 


tion 


M_^  • ABCj^  = Mj_  cos  0 


(3.2.19) 


and  is  equal  to  (Figure  3).  Solving  Equation  (3.2.19) 


for  r. 


Mj_  • ABCj^ 


(3.2.20) 


-Cjd  + . b^) 


/2.,j.,‘^2“‘^2.  /4._.  2.  2..  4..  2,.  2 


(3.2.21) 


The  preliminary  procedure  mentioned  in  Section  3 . 1 can 
now  be  carried  out.  Calculation  of  is  the  first  step. 


positive  sign  in  Equation  (3.2.15). 

The  angle  between  Mj_  and  ABCj_,  8»  satisfies  the  rela- 


At the  design  point,  5 and  0 are  the  stream  deflection  angle 
and  shock  angle  with  respect  to  the  free  stream  of  the  flow 
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“/vcr*-  . 


I 


I 


through  an  oblique  shock.  Thus  the  formula  for  the  Mach 

number  upstream  of  an  oblique  shock  can  be  used  to  find 

^des'  ^des  exist.  Equation  (3.2.22)  will  give 

a number  less  than  1. 

„2  ^ 10  (cot  9 -f  tan  5) 

”des  5 sin  2 0 - tan  <5(7  + 5 cos  20) 


The  second  step  is  to  determine  whether  or  not  the 
design  point  shock  is  weak  with  respect  to  the  free  stream. 
For  a given  Mach  number  the  shock  angle  can  range  from  the 
Mach  angle,  sin”^  to  90°,  depending  on  the  stream  deflec- 
tion angle.  At  the  maximum  stream  deflection  angle  the 
shock  angle  has  some  intermediate  value.  Below  this  are 
the  weak  shocks,  above  it  the  strong  ones  and  exactly  at 
this  value  the  shock  is  both  weak  and  strong.  Thus  the  de- 
sign point  shock  is  weak  with  respect  to  the  free  stream  if 
0 is  less  than  or  equal  to  the  shock  angle  which  would  give 
the  largest  possible  deflection  of  a streamline  at 

This  condition  is  expressed  in  Equation  (3.2.23). 


0 < sin 


-1 


7 

M^ 

des 

(3.2.23) 


■f 


ii! 
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The  third  step  is  to  determine  whether  or  not  $ > 't. 

— trans 


which  means  4'.  must  be  found.  Given  M,  ,6  and  0, 
trans  des 

4',, is  the  value  of  41  for  which  T - T,  goes  to  zero. 

trans  max  J- 

(5,0,$)  is  given  by  Equations  (3.2.21)  and  (3. 2. 4-6). 

^max  ® function  of  which  is  given  by 

Equations  (3.2.12),  (3.2.22),  (3.2.4)  and  (3.2.5).  The  rela- 
tion between  T and  M , is 

max  ± 


meoc 


= cot 


-1 


tan  ( ti);  p ) 
max 


6 M' 


5(M^sin^(u)  Ip  )-l) 
‘ max 


P .3.24) 


where  w|  is  the  shock  angle  that  would  occur  if  a 

mcuc 

streamline  were  deflected  through  the  angle  T and  is 

mcix 

given  by 


2 

sin  (wL  ) 
max 


3 Mf-5+  ^3(3M^  + 4 mJ  + 20) 
— 


(3.2.25) 


7 M 


The  equations  are  too  cxambersome  to  set  T - P equal  to 

niaX 

zero  and  solve  for  is  simpler,  especially  on 

the  computer,  to  search  for  the  value  of  $ at  which 
^max  ” zero.  Once  4*^^^^^^  is  found,  it  can  be 
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compared  to  i.  If  $ > , the  design  point  shock  is 

— trans 

weak  with  respect  to  M , and  calculation  of  , C_  and  L/D 

X Ij  D 

can  begin. 

3.3  Pressure  Ratio  Across  Lower  Surface  Shock 

An  attached  shock  will  not  be  possible  if  the  leading 
edge  is  subsonic  or  r |m  (M_,0,$)l  < r (6,0,<I>). 

(M^,0,'li)  is  given  by  Equations  (3.2.12),  (3.2.4)  and 
(3.2.5).  If  Mj^  < 1,  the  leading  edge  is  subsonic  and  calcu- 
lation of  the  aerodynamic  coefficients  for  the  current  values 

of  the  independent  variables  is  aborted.  Otherwise  T (M  ) 

max  j. 

is  calculated  using  Equations  (3.2.24)  and  (3.2.25)  and  is 
calculated  from  Equations  (3.2.21)  and  (3. 2. 4-6).  Again,  if 
< r further  calculation  is  terminated.  However,  if 

ntox  X 

^m2uc  — ® weak  shock  is  attached  to  the  leading  edge  and 

the  pressure  ratio  can  be  found  as  follows. 

No  simple  expression  for  the  pressure  ratio  in  terms  of 

and  Fj^  exists  so  it  is  necessary  to  calculate  the  shock 

angle,  u,  as  an  intermediate  step.  There  is  a cubic  equation 
2 

for  sin  w with  coefficients  that  depend  on  and  F^^ ; however, 
it  was  found  to  be  easier  to  se2u:ch  for  the  value  of  oj  be- 
tween sin  ^ (1/M  ) and  oj  1 - , which  gives  the  correct  value 

00  A 


of  F^  using  the  formula 


rj_  = cot  tan  0) 


5(M^sin^  oj-l) 


(3.3.1) 


Once  u is  known,  the  pressure  ratio  across  the  shock,  , 


is  given  by 


= 

p 

^ flO 


•7  «2  . 2 
7 M sxn  0)  -1 


(3.3.2) 


where  p^  and  p^  are  the  static  pressures  upstream  and  down- 
stream of  the  shock. 


3.4  Pressure  Ratio  Across  Upper  Surface  Expansion 


The  expansion  angle  normal  to  the  leading  edge, 

through  which  the  upper  surface  flow  turns  is  calculated  in 
the  same  manner  as  except  that  6',  c'  and  c^  replace 

<S,  C and  in  Equations  (3.2.3^,  (3.2.6)  and  (3.2.21). 
These  equations  become 


C = (l,c^,0) 


(3.4.1) 


= - tan  6 ' 


(3.4.2) 


= sin 


-c'  (1  + b^  + 


/c.2^1+(!^) 


/^b2+2b2b3+b3+b2+b3 


(3.4.3) 
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In  order  to  calculate  the  pressure  ratio  across  the 
upper  surface  expansion,  the  Mach  number  normal  to  the 

leading  edge  and  downstream  of  the  expansion,  must  be 

found.  To  find  the  Prandtl-Meyer  angle,  v,  through 
which  a streeun  would  turn  in  expanding  from  M-1  to 

is  needed  and  is  given  by 

v-2.4495  tan"^  ( .40825>il^-l)  -tan"^^M^-l  + . (3.4.4) 


is  found  using  the  following  equation  by  searching  for 
the  value  which  gives  the  correct  Prandtl-Meyer  angle,  v. 
v-2.4495  tan" ^ ( . 4 0 8 2 5^M^ ^ - 1) - tan” 1 (3.4.5) 


Once  M^^  is  )cnown,  the  pressure  ratio  can  be  determined 


with  the  formula 


5 - — 

u P« 


5 + M 


5 + M 


2 1 
u 

T 


-3.5 


(3.4.6) 


where  p^  and  p^  are  the  static  pressures  upstream  and  down- 
stream of  the  expansion. 
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3»5  Average  Skin  Friction  Coefficient 

The  calculation  of  four  quantities  will  be  described 
in  this  section.  They  are  the  average  laminar  and  turbulent 
skin  friction  coefficients  on  the  caret  wing's  lower  surface, 
and  and  the  average  laminar  and  turbulent  skin 


friction  coefficients  on  the  upper  surface,  and 

Reynolds  numbers  based  on  length  downstream  of  the 
lower  surface  shock  and  upper  surface  expansion  are  required 
for  the  calculation  of  the  skin  friction  coefficients.  A 
power  law  for  the  variation  of  viscosity  with  temperature 
of  the  form 


0.76 

^ = (|-)  , 300°R  < T < 900°R  (3.5.1) 

Wr  J-j- 


will  be  used  in  their  determination.  The  Reynolds  number 
ratio  across  the  lower  surface  shock  is  given  by 


where  ^2  length  of  the  internal  rib  and  the  sub- 

script I refers  to  a property  of  the  lower  surface  flow 
downstream  of  the  shock.  Expressions  for  the  ratios  in 
Equation  (3.5.2)  are 


J 


r / <t 

cos  o / tan  0 tan 


(3.5.3) 


/s 


p«  “ ® 


(3.5.4) 


1 


Mf(6 


1/2 


(3.5.5) 


0.76 


0.76 

6 +1  ' 


(3.5.6) 


Substituting  Equations  (3. 5. 3-6)  into  Equation  (3.5.2)  gives 


Re. 


- “ 

1 - 


5(5^1) 


M:(6ei+1) 


- A/  4. 

* 

[?7TqT5T 

0.76 


Re 


(3.5.7) 


cos  6 / tan  0 tan  *■ 
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The  Reynolds  number  ratio  across  the  upper  surface  expansion 
is  given  by 


21  _i 
/S 


(3.5.8) 


where  is  the  length  of  the  upper  surface  rib  and  the 

subscript  u refers  to  a property  of  the  upper  surface 
flow  downstream  of  the  expansion.  Expressions  for  the  ratios 
in  Equation  (3.5.8)  are 


cos  6'  /tan  0 tan 


(3.5.9) 


5 + M‘ 


5 + 


(3.5.10) 


5 + M' 


5 + m: 


(3.5.11) 


— = ( — ) 
u 'T  ' 

^u  u 


0.76 


5 + m" 


5 + m" 


0.76 


(3.5.12) 


where 


.-2/7 


(3.5.13) 


Substituting  Equations  (3.5.9-13)  into  Equation  (3.5.8)  gives 


_ Vcs+Mf)  - 5 

1 $ 

M cos  6 ' / tan  0 tan  ■=• 
00  ^ 


4.48 

7 


(3.5.14) 


The  Mach  number  of  the  lower  surface  flow  downstream 
of  the  shock  is  given  by 


M^(65j^+1)  - 5(?2-i) 


(3.5.15) 


The  four  average  skin  friction  coefficients  can  now  be 
calculated.  The  local,  laminar,  inccanpressible  skin  friction 
coefficient  is 


Cf^ (x) 


(3.5,16) 


where  x is  the  distance  from  the  leading  edge  in  the 
streamwise  direction.  The  average  over  the  lower  wetted 
surface  is  given  by 


**  \ A t 


dx  (x) 


2.656 


(3.5.17) 


where  h is  the  semi-span,  ^•(z)  = z (Figure  10)  and  the 

subscript  i indicates  incompressible  flow.  For  the  upper 
surface 


^^fiu^i 


2.656 


(3.5.18) 


Multiplying  Equations  (3.5.17)  and  (3.5.18)  by  the  ratio 
of  compressible-to-incOTipressible  laminar  skin  friction 
coefficients  gives 


2.656-.0718M^ 


> 0 < < 10 


(3.5.19) 


2.656-.0718M 


0 < < 10. 


(3.5.20) 


The  average  turbulent,  incompressible  skin  friction 


coefficient  between  the  leading  edge  and  a point  a distance 


I from  it  is 


C^.  (4)  » 5xl0^<Rej<10®. 

(log,o  Re,) 2-58  4 


(3.5.21) 


The  average  over  the  lower  wetted  surface  is 


CftC^-tz))  dz 


^/l 


-2.58 


log  Re,  + log 


dz  (3.5.22) 


where  l(z)  = z and  h is  the  semi-span  (Figure  10). 


Using  the  binomial  expansion 
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-2.58 


-3.58 


-2.58  (log  Re,  ) log  dz 


-2.58  5 -3 

.455  (log  Re,  ) + fe^dog  Re,  ) 

2 e 2 


.581 

L (3.5.: 


1 
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For  the  upper  surface 


'=ftu>i  = 


,455 


-2.58  , ..  -3.58" 

3 e 3 . 


(3.5.24) 


Multiplying  Equations  (3.5.23)  and  (3.5.24)  by  the  ratio  of 
compressible-to- incompressible  turbulent  skin  friction 
coefficients  gives 


( (1.0-0.1186  Mj^)(Cf^^)^,  0<Mj^<5.36 

tj?  “)  (3.5.25) 

( (0.5962-0.0430  Mj^)  (Cf^j^)^,  5.36<4-l^<10 


j (1.0-0.1186  M^^)  (Cf^^)^,  0<M^<5.36 
^ftu  ’ j 

( (0.5962-0.0430  M^) (C^^^)^,  5.36<M^<10 


(3.5.26) 


3.6  Base  Pressure  Ratio 

The  ratio  of  base  pressure  to  free  stream  static  pressure 
is  calculated  for  two  cases.  is  the  pressure  ratio  for  a 

turbulent  boundary  layer  approaching  the  trailing  edge  and 
is  the  pressure  ratio  for  a laminar  boundary  layer. 
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Base  pressure  ratio  is  correlated  to  the  ratio  of 
boundary  layer  thickness  at  the  trailing  edge  to  trailing 
edge  thickness.  This  parameter  depends  on  the  thickness 
ratio  and  Reynolds  number  based  on  chord  length.  The  thick- 
ness ratio,  T,  is  given  by 


T = (tan  6 - tan  6 ' ) cos 


(3.6.1) 


The  average  Reynolds  number  based  on  chord  length  over  the 
span  is 


Re. 


Re 


ave 


2 cos 


A 


tan  0 tan 


T 

I 


(3.6.2) 


For  a turbulent  boundary  layer  the  correlation 
parameter,  is  given  over  the  range  by 


T7T 


" <^«ave> 


(3.6.3) 


and  the  lines  fit  to  Chapman's  (11)  and  Goecke's  (12)  data 


are 
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/o. 5+0.0125 

(n^-0.5) , 

=1.5 

00 

jo. 35+0. 025 

(n^-0.5) , 

o 

. 

CM 

II 

8 

S 

jo.  2+0. 062 

(ti^-0.65)  , 

* 3.-1 

00 

\0. 14+0. 062 

(n^-0.22) , 

o 

• 

in 

li 

8 

S 

(3.6.4) 


At  intermediate  Mach  niimbers  linear  interpolation  between  the 
lines  is  used. 

For  a laminar  boundary  layer  the  correlation  parameter, 
is  given  over  the  range  0<^ri££0.2  by 


* 


1/2  • 


(3.6.5) 


Data  taken  from  Chapman's  (11)  curves  at  M^=1.5,  2.0  and  3.1 
is  presented  in  Table  1.  The  base  pressure  ratio  at  these 
three  Mach  numbers  is  found  by  parabolic  interpolation.  At 
intermediate  Mach  numbers  linear  interpolation  between  the 
curves  is  used. 

3.7  Lift  Coefficient,  Drag  Coefficient 
and  Lift-to-Drag  Ratio 

Formulas  for  the  aerodyneunic  coefficients  are  developed 
in  this  section.  These  expressions  depend  on  the  pressure 
ratios  and  skin  friction  coefficients  determined  in  the  pre- 
ceding sections,  the  ratio  of  free  stream  static  pressure  to 
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dynamic  pressure,  and  several  area  ratios  described 

in  this  section. 

Let  be  the  caret  wing  length  measured  from  the  apex 

to  the  plane  containing  the  base  (Figure  4) . Then  the  plan 
curea,  S,  is  given  by 

S * tan  0 tan  j . (3.7.1) 

The  area  of  the  projection  of  the  lower  surface  onto  a plane 
perpendiculeu:  to  the  free  stream  E^,  is  used  to  calculate  the 
contribution  of  the  lower  surface  press? .re  to  the  drag.  It 
is  given  by 

tan  0 tan  j tan  6 . (3.7.2) 

Similarly,  the  area  of  the  projection  of  the  upper  surface 
onto  a plane  perpendicular  to  the  free  stream,  E^,  is  used 
to  calculate  the  contribution  of  the  upper  surface  pressure 
to  the  drag.  Changing  6 to  <S'  in  Equation  (3.7.2)  gives 
the  following  expression  for  E^. 

E^  » tan  0 tan  j tan  6 ' 


(3.7.3) 
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The  base  area,  B,  is  just  the  difference  of  these  two  and  is 
given  by 

tan  0 tan  ^ (tan  6 - tan  6’)  . (3.7.4) 

The  area  of  a triangle  with  sides  of  length  a,  b and  c is 
given  by 

Area  » /s (s-a) (s-b)  (s-c)  (3.7.5) 

where  the  semi-perimeter,  s,  is 

s = j (a  + b + c)  . (3.7.6) 

The  wetted  eureas  of  the  lower  and  upper  surfaces,  Aj^  and  A^, 

are  each  composed  of  two  triangles.  For  the  lower  surface, 
the  two  triangles  have  sides  whose  lengths  are  given  by 


A S 

(3.7.7) 

a — 

cos  5 

b = 

Aj^y^an^Otan^ 

(tan0-tan5) ^ 

(3.7.8) 

c = 

ij^*^an^0tan^ 

tan^0+l  . 

(3.7.9) 

4 


!i 


The  wetted  area  of  the  lower  surface  is 


= 2 /s(s-a)  (s-b)  (s-c) 


(3.7.10) 


Changing  6 to  6'  in  Equations  (3.7.6-10)  gives  the 
wetted  area  of  the  upper  surface. 

Only  ratios  of  areas  appear  in  the  formulas  for 
C_ , C_  and  L/D,  so  Z.  drops  out. 

1j  U X 


These  ratios  are 


= tan  5 


-f  = tan  6- 


g'  “ tan  5 - tan  6 ' 


^ ^ 2 /s» (s'-a*) (s'-b*) (s'-c») 
® tan  0 tan  |- 


(3.7.11) 


(3.7.12) 


(3.7.13) 


(3.7.14) 


where 


'tan  0 tan  + (tan  0 - tan  6) 


/ 2 2~$  2 
'tan  0 tan^  j + tan  0 + 1 


= i (a'+b'+c') 


(3.7.15) 


(3.7.16) 


(3.7.17) 


(3.7.18) 
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Changing  <5  to  6'  in  Equations  (3.7.14-18)  gives  the  ratio 

V®- 

The  final  quantity  needed  in  the  formulas  for  the  aero- 
dynamic coefficients  is  the  ratio  of  free  stream  static 
pressure  to  dynamic  pressure.  It  is  given  by 


0 . 7 m' 


(3.7.19) 


Expressions  for  C_ , C and  L/D  can  now  be  derived.  The 

Xi  D 

lift,  L,  is  given  by 


L = (p^-p^j)S  - sin  6 - sin  5* 


(3.7.20) 


where  t . and  t are  the  boundary  layer  shear  stresses  on 
x»  u 

the  lower  and  upper  surfaces.  The  lift  coefficient  is  given 


q S ~ p q 

'loo  ~ 00  ^00 


The  drag,  D,  is  given  by 


T.  A„  T A 

Z Z . , u u . . 

— sin  5 - — — sin  6 


(C,)  ^ sin  S - (Cj)^  ^ sin  5' 


(3.7.21) 


i 


I 

I 
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D “ - Pu^u  " Pb®  ^ ■"  Vu  '^'  (3.7.22) 


where  p^  is  the  base  pressure. 


The  drag  coefficient  is  given  by 


'D  q«S 


n Poo 

- Su-T  - 5b  l>  =°=  5 


+ '=f>u  "I 


(3.7.23) 


The  lift-to-drag  ratio  is  given  by 


L/D  « • (3.7.24) 

» 

In  the  case  of  a laminar  boundary  layer,  | 

and  are  used  for  (Cj)^^,  (Cf)^  and  in  Equations 

(3.7.21),  (3.7.23)  and  (3.7.24).  For  a turbulent  boundary 
layer  case,  and  are  used  for  (Cj)j^,  (C^)^ 

and  For  the  case  of  zero  skin  friction,  (Cj)^^  and  (C^)^ 

are  set  to  zero.  For  zero  base  drag  is  set  equal  to  1.0. 

3.8  Computer  Program 

The  calculations  described  in  this  chapter  are  straight- 
forward but  quite  lengthy.  Therefore,  a computer  program  was 


written  for  the  calculation  of  the  aerodynamic  coefficients 
using  the  FORTRAN  IV  language.  Execution  of  this  program 
was  carried  out  at  M.I.T.'s  Information  Processing  Center 
on  their  IBM  System  370/168  computer.  Graphical  output 
was  produced  on  a CalComp  563  Drum  Plotter  using  a 
CalComp  905  Controller. 

The  computer  program  is  broken  up  into  subroutines 
along  the  lines  of  the  sectional  division  of  this  chapter. 
Subroutine  CARET  coordinates  the  basic  process  of  calculating 
C^  , C_  and  L/D  as  functions  of  6,6',0,'I>,M  and  Re  . Depend- 

^ Q 00  CO  * 

ing  on  the  test  matrix  desired,  a main  program  is  written 
to  feed  subroutine  CARET  with  a sequence  of  values  for  the 
independent  vauriables.  The  form  in  which  the  results  are 
recorded,  either  tabular  or  graphical,  is  also  controlled 
by  the  main  program.  A flow  chart  of  the  interaction  between 
the  main  program  and  subroutines  appears  in  Figure  11. 
Appendix  I contains  listings  of  the  main  programs  and  sub- 
routines and  a short  description  of  each. 
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CHAPTER  IV 

RESULTS 

4.1  Test  Matrix 

Testing  all  possible  combinations  of  the  six  inde- 
pendent variables  throughout  their  entire  ranges  would 
yield  an  enormous  and  unmanageable  amount  of  data.  There- 
fore, a test  matrix  was  devised  which  would  elucidate  most 
of  the  relationships  with  a reasonable  amount  of  computation. 
One  of  the  independent  variables,  5',  was  eliminated 
entirely  by  setting  it  equal  to  5-3°.  Essentially,  this 
restricts  the  study  to  wings  which  are  approximately  five 
per  cent  thick.  The  remaining  independent  variables  can  be 
broken  into  two  groups,  those  which  specify  the  caret  wing's 
geometry  and  incidence  eind  those  which  specify  the  flight 
conditions . The  test  matrix  was  divided  along  these  lines 
into  two  parts,  the  variable  geometry  test  and  the  fixed 
geometry  test. 

In  the  variable  geometry  test  the  flight  conditions, 
and  Re^,  were  held  fixed  while  the  geometric  p^u:ameters, 

5,  0 and  $,  were  allowed  to  vary.  This  test  is  composed  of 
three  parts  according  to  which  of  the  geometric  parameters 
is  being  examined.  When  the  dependence  of  the  aerodynamic 
coefficients  on  5 is  being  considered,  6 is  varied  from  3° 
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to  13°  in  steps  of  two  tenths  of  a degree.  When  0 is  being 
considered,  it  is  varied  from  4°  to  80°  in  steps  of  two  de- 
grees. When  i is  being  considered,  it  is  varied  from  24° 
to  178°  in  steps  of  two  degrees.  In  each  of  the  three  parts 
of  the  variable  geometry  test,  the  other  two  geometric 
parameters  are  allowed  to  take  on  all  combinations  of  the 
following  values:  6 =4°, 6°, 8°  and  10°;  0 = 30°,  50°  and  70°; 
and  * = 60°,  90°,  120°  and  150°. 

The  entire  variable  geometry  test  was  carried  out  at 
two  different  flight  conditions.  First  at  low  Mach  number 
and  high  Reynolds  number  with  a turbulent  boundary  layer: 

= 2.0  and  Re_  = 10^.  Second  at  high  Mach  number  and  low 
Reynolds  number  with  a leuninar  boundary  layer;  M^  » 3.1  and 

Re^  = 10^. 

00 

In  the  fixed  geometry  test  the  geometric  parameters, 

6,  0 and  $,  were  held  fixed  and  the  effects  of  Mach  number  I 

and  Reynolds  number  on  the  aerodynamic  coefficients  were 

examined.  This  test  is  composed  of  two  parts  according  to 

whether  the  skin  friction  and  base  drag  calculations  are 

based  on  a laminar  or  turbulent  boundary  layer.  For  the 

laminar  boundary  layer  case  the  Mach  number  is  varied  from 

1.5  to  3.1  by  tenths.  For  the  turbulent  boundary  layer  case 

the  Mach  number  is  varied  from  1.5  to  5.0  by  tenths.  Each 

part  of  the  fixed  geometry  test  is  repeated  for  Reynolds 

numbers  of  10^,  10®,  10^  and  10®. 
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The  particular  caret  wings  selected  for  the  fixed 
geometry  test  were  chosen  for  their  high  L/D  based  on  the 
results  of  the  variable  geometry  test.  For  the  laminar 
boundary  layer  case  the  caret  wing  <S  = 8.3°,  0 » 30.0°  and 
♦ ■ 120.0°  was  used.  For  the  turbulent  boundary  layer  case 
the  caret  wing  5 * 6.5°,  0 = 30.0°  and  <5  » 120.0°  was  used. 

Many  of  the  combinations  of  the  independent  variables 
which  occur  in  the  above  test  matrix  do  not  correspond  to 
caret  wings  that  belong  to  the  restricted  class  considered  in 
this  study.  They  may  be  disqualified  for  any  of  the  follow- 
ing reasons.  No  design  Mach  number  exists,  the  design  point 
shock  is  strong  with  respect  to  the  free  stream,  i is  less 


than 


, the  leading  edge  is  subsonic  or  T 

trans  max 


is  less 


than  . Such  cases  are  dropped  from  the  test  matrix  by 


the  computer  program  when  they  are  found. 

The  test  matrix  is  summarized  in  Table  2. 


4,2  Variable  Geometry  Test 


4.2.1  Dependence  on  6 


The  dependence  of  the  aerodynamic  coefficients  on 
6 is  illustrated  in  three  types  of  graphs.  Graphs  of  the 
first  type  indicate  the  effect  of  a second  geometric  parameter. 


This  is  done  with  a faunily  of  curves  for  four  different  values 
of  'I'  (60°,  90°,  120°  and  150°).  The  curves,  like  those  on 


all  graphs  described  in  this  chapter,  are  distinguished  from 
each  other  by  symbols  plotted  at  every  tenth  data  point. 

There  are  eighteen  graphs  of  this  type,  one  for  each  combina- 
tion of  the  three  aerodynamic  coefficients,  two  flight 
conditions  and  three  values  of  0 (30°,  50°  and  70°) . 

Graphs  of  C-  vs.  6,  vs.  6 and  L/D  vs.  5 appear  in 

D 

Figures  12-17,  24-29  and  36-41  respectively. 

Graphs  of  the  second  type  illustrate  the  contribution 
of  skin  friction  and  base  drag  to  the  overall  lift,  drag 
and  L/D.  The  three  curves  on  each  graph  correspond  to  in- 
clusion of  all  the  forces,  no  base  drag  and  no  base  drag  or 
skin  friction.  There  are  twelve  graphs  of  this  type,  one 
for  each  combination  of  the  three  aerodyncunic  coefficients, 
two  flight  conditions  and  two  values  of  0 (30°  and  50°) . 

The  value  of  $ is  120°  in  all  of  them.  Graphs  of  C_  vs.  6, 

Lf 

Cjj  vs.  6 and  L/D  vs.  6 appear  in  Figures  18-21,  30-33  and 
42-45. 

Graphs  of  the  third  type  illustrate  the  effect  of 
boundary  layer  type.  The  two  curves  on  each  graph  correspond 
to  lamineu:  and  turbulent  boundary  layers.  There  is  a graph 
for  each  combination  of  the  three  aerodynamic  coefficients 
and  two  values  of  0 (30°  and  50°) . All  six  of  these  graphs 

are  for  the  low  Mach  number,  high  Reynolds  number  flight  con- 
dition and  (f  equal  to  120°.  The  high  Mach  number,  low 


Reynolds  number  flight  condition  is  excluded  because  the 
Reynolds  number  is  outside  the  range  of  validity  of  the 
formula  for  the  turbulent  skin  friction  coefficient. 

Graphs  of  vs.  6,  Cj^  vs.  5 and  L/D  vs.  S appear  in 

Figures  22-23,  34-35  and  46-47. 

The  graphs  of  vs.  S appear  in  Figures  12-23.  As 

one  would  certainly  expect,  all  these  graphs  show  that  C_ 

Xi 

increases  with  5.  The  greater  the  value  of  6 , the 

greater  will  be  the  deflection  of  the  free  stream. 
Consequently,  the  pressure  ratio  across  the  lower  surface 
shock  will  increase  with  6 while  the  pressure  ratio  across 
the  upper  surface  expansion  decreases.  As  a result  the  lift 
coefficient  increases  with  6. 

Since  the  base  of  the  caret  wing  was  constructed 
perpendicular  to  the  free  stream,  the  base  pressure  ratio 
should  have  absolutely  no  effect  on  the  lift  coefficient. 

This  is  confirmed  by  Figures  18-21  in  which  the  curves  for 
all  forces  and  no  base  drag  lie  directly  on  top  of  each  other. 

Because  the  skin  friction  shear  stress  acts  tangent  to 
the  caret  wing's  surface  and  the  wing  is  at  fairly  low  angle 
of  attack,  skin  friction  should  have  only  a small  effect, 
which  increases  slightly  with  6 , on  the  lift  coefficient. 
Exactly  this  behavior  is  observed  in  Figures  18-21.  Since 
the  skin  friction  has  only  a small  effect  on  C_ , the  boundary 

li 
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layer  type  (Icuninar  or  turbulent)  should  not  make  a 
significant  difference  in  the  lift  coefficient.  This,  too, 
is  confirmed  by  Figures  22  and  23. 

Since  base  pressure  and  skin  friction  do  not  play  an 
important  role  in  the  behavior  of  the  only  forces  which 

make  a significant  contribution  are  those  due  to  the  upper 
and  lower  surface  pressures.  In  Figure  7 the  pressure  ratio 
across  the  lower  surface  shock  of  the  caret  wing  family 
6 = 5.0°,  0 = 34.0°  is  plotted  for  several  off  design  Mach 
numbers  as  a function  of  4.  When  < M^,  the  pressure 

ratio  increases  with  4 and  when  > M^,  the  press’^re 

ratio  decreases  with  increasing  4.  On  this  basis  one  might 
predict  that  will  increase  with  4 when  and 
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For  the  first  flight  condition,  M =2.0  and  Re  = 10^, 

the  lift  coefficient  ranges  from  a minimum  of  0.036  to  a 
meiximum  of  0.535.  For  the  second  flight  condition, 

= 3.1  and  Re^  = 10^,  ranges  from  a minimum  of  0.019 

to  a meucimum  of  0.313. 

The  graphs  of  Cj^  vs.  5 appear  in  Figures  24-35.  Many 
of  the  comments  made  about  the  dependence  of  C-  on  5 also 
apply  to  Cjj.  Again,  as  expected,  all  the  graphs  show  that 
Cjj  increases  with  6. 

Unlike  the  situation  for  the  base  pressure  and  skin 
friction  do  have  an  important  effect  on  Cj^.  For  the  flight 

7 

condition  = 2.0,  Re^  = 10  and  a turbulent  boundary  layer. 

Figures  30  and  31  show  that  the  base  drag  is  about  twice  the 
skin  friction  drag.  For  the  flight  condition  = 3.1, 

Re^  * 10^  and  a laminar  boundary  layer.  Figures  32  and  33 

show  that  the  skin  friction  drag  is  about  ten  times  the  base 
drag.  This  change  in  relative  importance  makes  sense  because 
base  drag  is  lower  for  a laminar  boundeury  layer  than  for  a 
turbulent  boundary  layer  and  skin  friction  drag  is  higher  for 
a lower  Reynolds  number.  For  both  flight  conditions. 

Figures  30-33  show  that  base  drag  and  skin  friction  drag  are 


the  principle  sources  of  drag  at  low  values  of  5 but  that 
pressure  drag  takes  over  at  higher  6.  Figures  34  and  35 
show  that  Cjj  is  higher  for  a turbulent  boundary  layer  than 
for  a laminar  boundary  layer  but  that  the  difference  is  small 
At  the  larger  values  of  6 , where  pressure  drag  is 
more  significant  than  base  drag  or  skin  friction  drag,  Cj^ 

displays  the  same  dependence  on  <t  as  discussed  for  C_  (see 

Li 

Figures  24-29) . At  lower  values  of  6 , the  pressure  drag 
becomes  insignificant  compared  to  the  other  sources  of  drag 
and  the  arguments  which  explain  the  dependence  of  on  4> 

can  no  longer  be  applied  to  Cj^.  As  $ increases,  there  is 

less  wetted  area  for  a given  plan  area  and  therefore,  the 
skin  friction  drag  goes  down  with  increasing  This  effect 

is  observed  in  Figures  25,  27  and  28  where  the  dependence  of 
Cjj  on  ♦ reverses  as  5 decreases  and  skin  friction  drag 

becomes  more  important  than  pressure  drag. 

For  the  first  flight  condition,  = 2.0  and  Re^  = 10^, 

the  drag  coefficient  ranges  from  a minimum  of  0.021  to  a 
maximum  of  0.133.  For  the  second  flight  condition,  * 3.1 
and  Re„  * 10^,  ranges  from  a minimum  of  0.022  to  a maximum 


of  0.093 
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The  graphs  of  L/D  vs,  <S  appear  in  Figures  36-47. 

Since  the  angle  6 is  closely  related  to  the  angle  of  attack, 
these  curves  would  be  expected  to  exhibit  a maximum.  Those 
curves  which  cover  a large  enough  range  of  5 have  the  ex- 
pected m^ucimum. 

Figures  42-45  show  the  contribution  of  skin  friction  and 
base  drag  to  I./D.  The  relative  importance  of  skin  friction 
and  base  drag  is  the  same  here  as  described  for  the  drag 
coefficient.  For  M » 2.0  and  Re  =*  10  , base  drag  is 
significant  compared  to  skin  friction.  For  = 3.1  and 
Re^  = 10^,  base  drag  is  insignificant  compared  to  skin  fric- 
tion. Figures  46  and  47  show  that  the  higher  drag  associated 
with  a turbulent  boundary  layer  leads  to  a lower  L/D  than 
for  a laminar  boundary  layer. 

Figures  37-41,  in  which  show  a predcsninant 

tendency  for  L/D  to  increase  with  $.  However,  there  are 
sane  exceptions  in  Figures  39  and  40.  Figure  36,  in  which 
Mdes  ^ shows  the  opposite  trend.  In  this  figure  L/D  de- 
creases with  increasing  4. 

For  the  first  flight  condition,  M =2,0  and  Re  = 10^, 

the  meucimum  L/D  occurs  at  6 » 6.4°  on  the  curve  $ = 120°  and 
is  equal  to  5.679,  For  the  second  flight  condition,  = 3.1 

and  Re  = 10^,  the  maximum  L/D  occurs  at  6 * 8.4°  on  the 
curve  ♦ » 120°  and  is  equal  to  4.057. 
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4.2.2  Dependence  on  0 

The  dependence  of  the  aerodynamic  coefficients 
on  0 is  illustrated  in  Figures  48-65.  These  figures 
also  indicate  the  dependence  on  6 with  a family  of 
curves  corresponding  to  four  values  of  that  parameter 
(4°,  6°,  8°  and  10°).  There  is  a graph  for  every  combina- 
tion of  the  three  aerodynamic  coefficients,  two  flight 
conditions  and  three  values  of  $ (90°,  120°  and  150°) . 

The  design  point  has  been  marked  on  the  four  curves  of 
Figure  48  and  occurs  at  the  same  values  of  0 in 
Figures  49-65. 

In  Section  4.2.1  it  was  established  that  the  lift 
coefficient  is  primarily  influenced  by  upper  and  lower 
surface  pressure  rather  than  skin  friction  or  base 
pressure.  Examining  the  expression  for  Cj^,  Equation 
(3.7.21),  these  pressures  enter  in  terms  of  the  pressure 
ratios  across  the  lower  surface  shock  and  upper  surface 
expansion  in  such  a way  that  increases  with  increasing 
and/or  decreasing  C . The  pressure  ratios  are 
functions  of  Mj_,  and  , which  in  turn  depend  on  0. 

Sample  calculations  within  the  ranges  of  6,6'  and 
4 of  interest  show  that  F^^  and  increase  with  decreasing 
0.  Since  increases  with  Fj_  and  decreases  with 
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increasing  , it  follows  that  will  tend  to  increase 

with  decreasing  0 due  to  changes  in  the  stream  deflection 
angles. 

Combining  Equations  (3.2.4),  (3.2.5)  and  (3.2.12) 
gives  an  expression  for  Mj_  as  a function  of  0.  The 
derivative  of  Mj_  with  respect  to  0 is  always  positive, 
so  increases  with  0.  However,  the  pressure  ratios 

across  the  shock  and  expansion  may  increase  or  decrease 
with  increasing  Mj_  depending  on  the  values  of  Mj_,  rj_ 
and  T2 . Thus  the  dependence  of  on  0 due  to  changes 

in  the  normal  component  of  the  Mach  number  will  vary  with 
the  situation. 

In  light  of  these  results,  the  overall  dependence  of 
on  0 has  not  been  anticipated.  Looking  at  the 
graphs  of  vs.  0 which  appear  in  Figures  48-53,  it 

can  be  seen  that  C.  increases  with  decreasing  0 through- 
out  the  range  of  the  present  test  matrix.  The  graphs  are 
also  consistent  with  the  previous  finding  that  C_  increases 

L 
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with  6.  For  the  flight  condition  = 2.0,  Re  = 10  , the 
lift  coefficient  ranges  from  a minimum  of  0.077  to  a maximum 
of  0.484.  For  the  flight  condition  = 3.1,  Re^  = 10^, 
the  lift  coefficient  ranges  from  0.045  to  0.263. 
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The  remarks  made  about  the  dependence  of  C_  on  0 

Li 

due  to  upper  and  lower  surface  pressure  also  apply  to  the 
drag  coefficient.  The  changes  in  and  C,,  which 

cause  to  decrease  with  increasing  0 will  also  cause 

the  pressure  drag  to  decrease  with  increasing  0. 

Besides  pressure  drag,  skin  friction  and  base  drag 
also  influence  the  drag  coefficient.  The  ratio  of  base 
eurea  to  plan  area  given  by  Equation  (3.7.13)  does  not 
depend  on  0.  However,  the  Reynolds  number  based  on 
length  decreases  with  increasing  0 leading  to  a higher 
base  pressure  ratio  and  lower  base  drag.  Thus  base  drag 
reinforces  the  effect  of  pressure  drag,  driving  the  drag 
coefficient  down  when  0 increases. 

Equations  (3.7.14-18)  show  that  the  ratio  of  wetted 
area  to  plan  area  increases  with  0.  This  means  that  skin 
friction  drag  and  its  importance  relative  to  the  other 
sources  of  drag  increase  with  0.  At  large  enough  values 
of  0 skin  friction  drag  might  predominate  and  reverse 
the  dependence  of  on  0. 

This  reversal  is  observed  clearly  at  the  high  Mach 
number,  low  Reynolds  number  flight  condition  in  Figures  57-59. 
The  curves  of  Cj^  vs.  0 have  a negative  slope  at  low  values 
of  0 and  positive  slope  at  high  values  of  0.  It  was 
shown  in  Section  4.2.1  that  the  magnitude  of  the  pressure 
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drag  relative  to  skin  friction  drag  increases  with  5.  This 
accounts  for  the  shift  of  the  minimum  in  these  curves  to  the 
right  with  increasing  6. 

It  was  also  shown  in  Section  4.2.1  that  base  drag  is 
greater  than  skin  friction  drag  at  the  low  Mach  number, 
high  Reynolds  number  flight  condition.  Most  of  the  curves 
corresponding  to  this  flight  condition,  Figures  54-56,  do 
not  have  a minimum.  Except  for  a few  cases  at  low  5 and 
high  0,  the  drag  coefficient  decreases  with  increasing  0. 

The  graphs  for  both  flight  conditions  are  consistent 
with  the  previous  finding  that  increases  with  5. 

For  M_^  ■ 2.0  and  Re^  ■=  10^  the  drag  coefficient  ranges 
from  a minimum  of  0.024  to  a maximum  of  0.097.  For 
M^  ■ 3.1  and  Re^  ■ 10^  the  drag  coefficient  ranges  from 
0.021  to  0.066. 

J 

The  dependence  of  the  lift-to-drag  ratio  on  0 is 
shown  in  Figures  60-65.  Due  to  the  effect  of  skin  friction, 

C_  does  not  fall  off  as  rapidly  as  C.  when  0 is 

U ii 


increased.  As  a result,  the  lift-to-drag  ratio  decreases 
with  increasing  0.  For  the  M^  » 2.0,  Re_j^  • 10^  flight 
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6 = 6°,  0 = 16°  and  f = 120°  and  has  the  value  4.973. 

Figures  60-65  are  consistent  with  the  previous  finding 
that  the  lift-to-drag  ratio  goes  through  a meiximum  as  5 
increases . 

4.2.3  Dependence  on  i 

The  dependence  of  the  aerodynamic  coefficients 
on  $ is  illustrated  in  Figures  66-83.  These  figures 
also  indicate  the  dependence  on  5 with  a family  of  curves 
corresponding  to  four  values  of  that  parameter  (4°,  6°,  8° 
and  10°) . There  is  a graph  for  every  combination  of  the 
three  aerodynamic  coefficients,  two  flight  conditions 
and  three  values  of  0 (30°,  50°  and  70°) . 

Figures  12-17,  24-29  and  36-41  discussed  in  Section  4.2.1 
have  already  indicated  the  dependence  of  C-  , C^.  and  L/D  on  $ . 

Xj  u 

Essentially  the  same  dependence  is  observed  in  Figures  66-83. 

It  was  established  in  Section  4.2.1  that  the  lift  coefficient 
increases  with  4>  when  and  decreases  with  in-  '■* 

creasing  $ when  ^ I*'  Figures  67-71  is  less 

than  and  C_  is  observed  to  increase  with  4'  in  most 

cases.  A few  exceptions  are  present  in  which  C_  begins 

Xj 

to  decrease  as  4>  approaches  180°.  In  Figure  66  is 

greater  than  and  decreases  with  increasing  $.  j 
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At  = 2.0  and  Re^  = 10^  the  lift  coefficient  ranges 
from  0.051  to  0.484.  At  M =3.1  and  Re  = 10^  C.. 

ranges  from  0.026  to  0.225. 

The  drag  coefficient  shows  the  same  dependence  on  4* 
as  C_  except  at  low  enough  values  of  5 and  i where 
the  effect  of  skin  friction  can  reverse  the  behavior 
(Figures  72-77).  Cj^  ranges  from  0.024  to  0.097  for  the 
first  flight  condition  and  from  0.023  to  0.091  for  the 
second  flight  condition. 

The  increase  in  lift-to-drag  ratio  with  4>  when 
Mdeg  observed  in  Section  4.2.1  is  confirmed  by 

Figures  79-83.  Exceptions  to  this  trend  which  were 
noted  in  that  section  are  seen  here  as  a sharp  decrease 
in  L/D  as  i approaches  180°.  Figure  78  is  consistent 
with  the  previously  observed  decrease  in  L/D  with 
increasing  $ when  > M^.  The  meucimum  L/D  for 

the  first  flight  condition  occurs  at  6 = 6°,  0 = 30° 
and  4>  = 100°  and  has  the  value  6.006.  The  maximum  L/D 
for  the  second  flight  condition  occurs  at  6 = 8°,  0 = 30° 
and  ♦ = 124°  and  has  the  value  4.050. 

The  dependence  of  C- , C_  and  L/D  on  6 in  Figures 

ij  u 

66-83  is  the  Scune  as  described  in  Sections  4.2.1  and  4.4.4. 
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4.3  Fixed  Geometry  Test 

The  dependence  of  the  aerodynamic  coefficients  on 
Mach  number  and  Reynolds  number  is  illustrated  in 
Figures  84-89.  Mach  number  is  the  abscissa  and  Reynolds 
number  a parameter  in  these  figures.  There  is  a graph 
for  every  combination  of  the  three  aerodynamic  coefficients 
and  two  boundary  layer  types.  The  caret  wing  used  in  the 
turbulent  boundary  layer  case  is  on  design  at  M^  = 2.38 
and  the  one  used  in  the  laminar  boundary  layer  case  is  on 
design  at  M^  = 2.53. 

In  linearized  thin  airfoil  theory  (18)  the  dependence 

of  the  pressure  coefficient  on  free  stream  Mach  number  is 

2 -1/2 

proportional  to  (M^  - 1)  . Since  the  deflection  of  the 

streamlines  in  the  caret  wing  flow  field  is  small,  the 
dependence  of  the  lift  and  drag  coefficients  on  Mach 
number  should  be  fairly  similar  to  this.  As  a further 
indication  of  the  influence  of  Mach  number,  consider  the 
terms  in  the  expressions  for  and  C_  , Equations 

(3.7.21)  and  (3.7.23).  The  difference  between  the  pressure 
ratios  across  the  shock  and  expansion,  multiplied  by  the 
ratio  of  free  stream  static-to-dynamic  pressure  is  propor- 
tional to  the  difference  in  pressure  coefficient  between 
the  upper  and  lower  surfaces.  For  small  stream  deflection 
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angles  this  term  is  proportional  to  (M^  - 1)  , 

corresponding  to  thin  airfoil  theory.  As  described  in 
Sections  2.4  and  3.5,  a linear  approximation  to  the  de- 
crease in  skin  friction  coefficient  with  increasing  Mach 
number  is  used.  The  contribution  of  base  drag  to  Cj^ 
is  proportional  to  (1  - Cjj)Poo/<3«*  The  factor  (1  - 

increases  slowly  with  but  the  ratio  of  static-to- 

_2 

dynamic  pressure  falls  off  with  . Overall,  this  term 
decreases  as  increases. 

These  remarks  are  consistent  with  the  dependence  of 
and  Cp  on  shown  in  Figures  84-87.  The  curves 

of  L/D  VS.  for  the  turbulent  boundary  layer  case. 

Figure  88,  have  a minimum  near  = 3.2.  The  highest 
L/D  on  each  curve  occurs  at  the  low  end  of  the  Mach  number 
range.  For  the  laminar  boundary  layer  case.  Figure  89, 
the  maximum  L/D  on  each  curve  also  occurs  at  the  low  end 
of  the  Mach  number  range.  Two  of  the  curves  decrease 
throughout  the  range  and  two  have  a minimum. 

Reynolds  ntunber  influences  the  aerodynamic  coefficients 
through  the  skin  friction  coefficients  and  base  pressure 
ratio.  The  skin  friction  drag  decreases  and  the  base  drag 
increases  as  Re_  increases.  In  most  cases  the  skin  friction 
effect  is  the  primary  one. 
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Changes  in  skin  friction  have  little  effect  on  the  lift 
coefficient  and  changes  in  base  drag  have  no  effect  so 
should  be  nearly  independent  of  Re^.  The  curves  for 
different  Reynolds  numbers  in  Figures  84  and  85  confirm 
this  by  lying  virtually  on  top  of  each  other. 

For  a turbulent  boundary  layer.  Figure  86  shows  that 
the  drag  coefficient  decreases  as  Re  increases.  This 
is  also  true  in  the  laminar  boundary  layer  case  for  Mach 
numbers  above  2.3  (Figure  87).  For  lower  Mach  numbers 
there  is  a narrow  range  of  Reynolds  numbers,  between  10^ 
and  10^,  in  which  the  effect  of  base  drag  dominates  and 
the  dependence  of  Cj^  on  Re„  reverses.  This  coincides 
with  the  steepest  portion  of  the  laminar  base  pressure 
curves  in  Figure  9. 

The  dependence  of  the  lift-to-drag  ratio  on  Re^ 
is  just  the  inverse  of  the  drag  coefficient's  behavior. 

For  a turbulent  boundary  layer  L/D  increases  with 
Re^  (Figure  88).  This  is  also  true  in  the  laminar 
boundary  layer  case  above  M^  = 2.3  (Figure  89).  At  lower 
Mach  numbers,  the  lift-to-drag  ratio  decreases  over  some 
range  of  Reynolds  numbers  between  10®  and  10^. 
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CHAPTER  V 

CONCLUSIONS  and  RECOMMENDATIONS 

A clear  picture  of  the  way  the  lift  coefficient,  drag 
coefficient  and  lift-to-drag  ratio  vary  with  caret  wing 
geometry  and  flight  condition  has  emerged.  The  influence 
of  6 , 0 , ♦ , and  Re^  on  each  of  the  aerodynamic  co- 
efficients is  summarized  below. 

The  lift  coefficient  increases  with  5 under  all 

conditions  and  with  when  M.  < M . It  decreases  as 

des  “ 

0 or  M increase  and,  when  M,  > M , as  $ increases. 

00  ' aes  ” 

Reynolds  number  has  no  appreciable  effect  on  C_ . 

The  drag  coefficient  also  increases  with  5 under 
all  conditions.  Drag  due  to  upper  and  lower  surface 
pressure  decreases  as  0 increases.  However,  skin  fric- 
tion drag  increases  with  0 and  can  be  the  dominant 
factor  at  some  flight  conditions.  Depending  on  whether 
^des  greater  than  or  less  than  M^,  the  drag  coefficient 
will  decrease  or  increase,  respectively,  as  increases. 
This  dependence  on  arises  through  the  influence  of  upper 
and  lower  surface  pressure.  At  low  values  of  5 and  <1> 
the  effects  of  skin  friction  may  override  this  dependence 
and,  in  such  cases,  will  decrease  as  $ Increases  regard- 
less of  the  relative  magnitudes  of  and  M^.  Finally,  the 


r 
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drag  coefficient  decreases  as  M or  Re  increase.  The 
exception  to  this  dependence  on  Re^  noted  in  Section  4.3 
is  a consequence  of  the  rapid  rise  in  lamincu:  base  drag 
over  a narrow  Reynolds  number  range.  This  result  must  be 
viewed  with  suspicion  because  the  two-dimensional  base 
pressure  data  gives,  at  best,  only  a rough  indication  of 
the  base  drag. 

The  lift-to-drag  ratio  goes  through  a mciximum  as  6 

is  Vcuried  and,  when  M, < there  is  a meucimum  in  its 

des  <*> 

*1 

dependence  on  When  ^ decreases  as  i 

increases.  The  lift-to-drag  ratio  also  decreases  with 
increasing  0.  For  cases  where  the  Mach  number  range 
extended  up  to  5.0,  a minimum  was  found  in  the  dependence 
of  L/D  on  M^.  Finally,  the  lift-to-drag  ratio  increases 
with  The  remarks  made  about  the  exception  to  the 

dependence  of  on  Re^  also  apply  here. 

The  best  lift-to-drag  ratio  found  in  the  variable 
geometry  test  was  6.483  for  the  caret  wing  5 = 6°,  0 = 14° 
and  $ * 150°  and  the  flight  condition  M = 2.0,  Re  = 10^. 

The  corresponding  lift  coefficient  and  design  Mach  number 
are  0.285  and  5.989.  Taking  the  best  L/D  from  each  of  the 
graphs  of  L/D  vs.  5,  0 and  *,  two  other  cases  were  found 
in  which  the  lift-to-drag  ratio  exceeded  6.0  and  many  were 
found  between  5.0  and  6.0. 


The  best  L/D  found  in  the  fixed  geometry  test  was 
6.391  for  the  caret  wing  5 = 6.5°,  0 = 30°  and  $ = 120° 

g 

and  flight  condition  = 1.8,  Re^  = 10  . The  correspond- 
ing lift  coefficient  and  design  Mach  number  are  0.324  and 
2.384. 

These  examples  of  high  L/D  are  only  the  best  cases 
from  the  test  matrix  considered  in  this  study.  Some  im- 
provement can  be  expected  by  searching  in  the  neighborhood 
of  these  cases,  using  the  behavior  of  L/D  described  above 
as  a guide.  Additional  improvement  in  L/D  could  be 
achieved  by  adding  scane  type  of  afterbody  and  thereby 
reducing  the  base  drag.  The  curves  for  zero  base  drag 
in  Figures  42-45  give  an  indication  of  the  room  for  improve 
ment.  In  Figure  42,  for  example,  the  L/D  goes  from  5.679 
to  8.945  when  base  drag  is  eliminated. 

The  performance  of  caret  wings  is  compared  to  that  of 

delta  wings  in  Figure  90.  Wind  tunnel  tests  of  8 percent 

thick  delta  wings  with  double  wedge  profiles  at  = 1.92 

and  Reynolds  number  based  on  two-thirds  the  root  chord,  Re^ 

6 6 

covering  the  range  0.57  x 10  to  1.25  x 10  were  made  by 
Love  (19).  A thinner  series  of  delta  wings,  1.3  to  1.8 
percent  thick,  constructed  from  flat  plate  with  sharpened 
leading  and  trailing  edges  are  also  included  in  his  tests. 
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The  boundary  layer  over  most  of  the  wetted  surface  was 
shown  to  be  turbulent  by  a liquid-film  method.  Love 
measured  the  aerodynamic  coefficients  through  an  angle 
of  attack  range  of  ± 6°.  The  meucimum  L/D  of  these  delta 
wings  is  plotted  as  a function  of  apex  half-angle  in 
Figure  90. 

The  free  stream  conditions  in  the  first  part  of  the 

7 

variable  geometry  test,  = 2.0  and  Re^  = 10  , are 
nearly  the  same  as  the  conditions  in  Love's  experiment. 
Ceuret  wings  with  high  L/D  from  this  part  of  the  study  are 
plotted  alongside  the  delta  wing  data  in  Figure  90  using 
the  apex  half-angle  of  the  planform  for  the  abscissa. 

The  Mach  number  and  average  Reynolds  ninnber  of  Love's 
experiment  were  duplicated  exactly  at  one  point  of  the 
test  matrix  in  the  fixed  geometry  test  and  this  case  is 
also  included  in  the  figure.  The  caret  wing  data  does  not 
represent  the  maximum  L/D  for  any  given  planform,  only  the 
best  to  occur  in  the  test  matrix.  Geometric  properties, 
and  Re-  are  listed  in  tabular  form  on  the  figure. 

The  thickness  ratio  of  the  caret  wings  is  .052,  about 
midway  between  the  two  groups  of  delta  wings.  As  shown  in 
Figure  90,  the  lift-to-drag  ratios  of  the  c2Lret  wings  also 
lie  between  the  two  groups  of  delta  wings,  indicating  that 
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caret  and  delta  wings  have  comparable  lifting  efficiencies. 
However,  the  delta  wings  considered  here  have  sharp  trail- 
ing edges  while  the  caret  wings  have  blunt  bases.  If  the 
addition  of  an  afterbody  to  the  caret  wings  could  sig- 
nificantly reduce  their  base  drag,  the  resulting  lift-to- 
drag  ratios  would  be  superior  to  those  for  delta  wings. 

One  objective  of  further  reseeurch  should  be  the 
determination  of  the  accuracy  and  range  of  applicability 
of  the  method  of  calculation  presented  in  this  thesis. 

This  could  be  accomplished  through  experimental  measure- 
ments or  calculations  based  on  a more  sophisticated  model 
of  the  caret  wing  flow  field. 

The  work  of  Ziph  (7)  mentioned  in  the  Introduction 
could  be  used  as  a steurting  point  in  the  latter  approach. 
Ziph  treated  the  case  of  a reentrant  pyramid  with  surface 
curvature  in  the  transverse  plane.  His  numerical  solution 
of  the  equations  of  conical  flow  is  restricted  to  flow 
fields  in  which  the  shock  is  attached  to  the  apex  but  de- 
tached from  the  leading  edges.  This  restriction  is  the 
result  of  a transformation  of  the  region  of  computation 
into  a square.  The  computational  region  is  bounded  by  the 
shock,  the  body  and  two  symmetry  lines  extending  radially 
from  the  leading  edge  and  internal  rib  to  the  shock.  If 
the  shock  were  attached,  the  length  of  one  side  of  the 
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computational  region  would  go  to  zero  and  the  transformation 
would  have  a singularity. 

As  shown  in  Chapter  II/  the  leading  edge  shock  can  be 
calculated  correctly  by  treating  the  caret  wing  as  a swept 
wedge.  One  way  of  extending  Ziph's  method  to  the  attached 
shock  case  would  be  to  modify  the  computational  region  so 
that  the  part  near  the  leading  edge  is  no  longer  included. 
The  leading  edge  shock  of  Chapter  II  would  then  be  used 
in  this  region  auid  would  give  the  boundary  conditions 
necessary  for  the  new  computational  region.  Neglecting 
boundary  layer  effects,  this  procedure  would  yield  the 
correct  shock  shape  and  pressure  distribution  over  the 
caret  wing's  lower  surface.  The  error  incurred  by 
approximating  the  average  pressure  over  the  lower  surface 
with  the  pressure  downstream  of  the  leading  edge  shock 
could  then  be  assessed.  If  found  to  be  unacceptably 
large,  this  approximation  could  be  abandoned  in  favor  of 
Ziph's  lengthier,  but  more  accurate  method. 

Another  objective  of  further  research  would  be  the 
inclusion  of  caret  wings  whose  design  point  shocks  are 
strong  with  respect  to  flow  normal  to  the  leading  edge 

) in  the  search  for  high  lift-to-drag  ratios. 

trans 


Jl 
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If  these  "strong"  shocks  actually  occur,  the  high  lower  : ; 

surface  pressure  that  could  be  generated  at  low  angle  of 
attack  would  yield  lift-to-drag  ratios  much  higher  than 
those  found  in  this  study.  An  experimental  investiga- 
tion of  this  possibility  would  be  very  interesting. 
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APPENDIX  I 

COMPUTER  PROGRAM  LISTING 

A short  description  of  the  subroutines  precedes  the 
listing. 

Subroutine  HEAD  prints  a heading  for  data  printed  by 
subroutine  PRINT.  There  are  no  arguments. 

Subroutine  CARET  calculates  design  point  shock 

type  with  respect  to  flow  normal  to  the  leading  edge, 

^u'  ^bil'  ^bt'  ^L'  ^D  given  5,  6’,  0,  $,  M^  and  Re^. 

Subroutine  PRINT  prints  data  under  the  heading  pro- 
duced by  subroutine  head . Input  arguments  are  6 , 6 ' , 8 , 
ir  Re„,  M^^g,  shock  type,  C^,  C^,  C^^ 

and  L/D. 

Subroutine  LOWER  calculates  M^gg/  shock  type, 

^ftA  5iven  5,  0,  $,  M^  and  Re^. 

Subroutine  DESIGN  calculates  and  given  6 

and  0 . 

Subroutine  DG  calculates  T T.  , - F , and  M. 

max  J-  max  i. 

given  5,  0,  <f  and  M^.  It  also  calculates  1^'  when  given  5' 
instead  of  6 . 

Subroutine  SHOCKW  calculates  the  pressure  ratio  across 
a weak  shock  given  the  upstream  Mach  number  and  stream  de- 
flection angle. 
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Subroutine  SHOCKS  calculates  the  pressure  ratio  across 
a strong  shock  given  the  upstream  Mach  number  and  stream 
deflection  angle. 

Subroutine  SFL  calculates  and  given  5,  0, 

Re*  and 

Subroutine  UPPER  calculates  *“f£u  ^ftu 

6 * , 0 / i,  M and  Re  . 

00  00 

Subroutine  PM  calculates  the  pressure  ratio  across  a 
Prandtl-Meyer  expansion  given  the  upstream  Mach  number  and 
angle  through  which  the  stream  turns. 

Subroutine  SFU  calculates  and  given  S',  0, 

i,  K,  Re„  and  5,  . 

Subroutine  BASE  calculates  and  given  S,  S', 

Q,  M_j^  and  Re^. 

Subroutine  COEFF  calculates  , C_  and  L/D  given  5 , 6 ' , 

li  u 

Qf  and 

Subroutine  PICTR  draws  labeled  axes,  title  and  sub- 
title and  plots  up  to  four  curves  identified  by  symbols. 
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SIIDROtiTINE  DFSIGN  (DRLT  A , T HETA  , XHDES  , PHIB) 

DELTA  = ANGLE  BETWEEN  FFEE  STREAM  6 INTERNAL  RIB,  IN  DEGREES 

THETA  = ANGLE  BETWEEN  FREE  STREAM  6 PLANE  SPANNING  L.E.'S,  IN  DEGREES 

XMDES  = DESIGN  MACH  NC. 

PHIB  = VALUE  OF  PHI  WHICH  DIVIDES  WEAK  6 STRONG  SHOCK  CASES  (ON  DESIGN) 


SUBflOUTTNH  DG ( DELT A , TH ETA , PHI , XM , DM, G AMM A , DMG , XH) 

DELTA  = ANGLE  BETWEEN  FREE  STREAM  R INTERNAL  RIB,  IN  DEGREES 

THFTA  = ANGLE  BETWEEN  FREE  STREAM  S PLANE  SPANNING  L.E.'S,  IN  DEGBEFS 

PHI  = ANGLE  BETWEEN  FREE  STREAM  SURFACES  CONTAINING  LEADING  FDGBS,  CEGREFS 
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SUBBCUTINE  SHOCK W ( X N , GKMH A , f R ATIO) 

XN  = HACK  NUHEFR 

GAHMA  = DEFLECTION  ANGLE,  IN  DEGREES 
PRATIO  = PRESSOHE  RATIO  ACROSS  SHOCK 
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SUBROUTINE  SHOCKS  (X N , 0 AN H A , ER AT 10) 

XN  = NACH  NUNPER 

r.ABHA  = DEFLECTION  ANGIE,  IN  DEGREES 
PPATIO  = PRESSURE  PATIO  ACROSS  SHOCK 
PPATIO  SET  = C.  IF  XN  .LE.  1.  OR  SHOCK  IS  DETACHED 
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Table  1 


LAMINAR  BASE 

PRESSURE  DATA 

^bJl 

n . 

M =1.5 

M =2.0 

M » 3 

i 

00 

00 

00 

0.00 

0.494 

0.338 

0.160 

0.01 

0.444 

0.357 

0.257 

0.02 

0.557 

0.570 

0.376 

0.03 

0.647 

0.669 

0.472 

0.04 

0.704 

0.716 

0.554 

0.05 

0.738 

0.742 

0.631 

0.06 

0.764 

0.764 

0.697 

0.07 

0.782 

0.779 

0.726 

0.08 

0.800 

0.789 

0.09 

0.809 

0.795 

0.10 

0.822 

0.800 

0.11 

0.825 

0.800 

0.12 

0.831 

0.806 

0.13 

0.838 

0.809 

0.14 

0.844 

0.813 

0.15 

0.844 

0.813 

0.16 

0.844 

0.816 

0.17 

0.847 

0.819 

0.18 

0.847 

0.819 

0.19 

0.847 

0.819 

0.20 

0.847 

0.819 

Table 


Figure  1.  Construction  of  Caret  Wing  from  Knovm  Flow  Field 


Figure  2.  Classification  of  Caret  Wings 


Figure  4.  Family  of  Caret  Wings  Derived  from  a Single  Wedge  Flow 
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DRAG  RRTIO  VS.  THETfl  Figure  62. 
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RRTIO  VS.  PHI  Figure  82. 
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Figure  87. 
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